Acousto-Optical Coherence Tomography with a digital holographic 

detection scheme 
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Acousto-Optical Coherence Tomography (AOCT) consists in using random phase jumps on ultrasound and 
light to achieve a millimeter resolution when imaging thick scattering media. We combined this technique with 
heterodyne off-axis digital holography. Two-dimensional images of absorbing objects embedded in scattering 
phantoms are obtained with a good signal to noise ratio. We study the impact of the phase modulation 
characteristics on the amplitude of the acousto-optic signal and on the contrast and apparent size of the 
absorbing inclusion. © 2012 Optical Society of America 

OCIS codes: 110.0113, 110.7170, 170.3880, 170.3660, 090.2880, 090.1995. 
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Medical imaging is a contrast issue. Different types of 
waves (X-rays, ultrasound, etc.) are used depending on 
which organ is observed. Imaging with light, e.g. for 
breast cancer screening, raises the problem of detect- 
ing a millimeter-sized absorbing object in a several- 
centimeters thick scattering medium. When using light 
only, like in Diffuse Optical Tomography [1], the resolu- 
tion of a breast tissues image is usually around 10 mm 
so that doctors can hardly detect an emerging tumor. 
Even if another imaging technique is used for screening, 
like radiography, the optical information remains pre- 
cious since it supplements the knowledge of the tissues 
with rich physiological details [2] . The idea to use acous- 
tic waves to achieve a millimeter resolution on optical 
information [3] gave rise to Ultrasound modulated Opti- 
cal Tomography (UOT). 

UOT is based on the acousto-optic (AO) effect which 
enables to tag scattered light by shifting its frequency lot, 
of the acoustic frequency lotjs [4]. The detection of this 
weak signal at lot, ± lotjs provides spatially resolved op- 
tical information by scanning the acoustic source within 
the medium. The size of the tagging zone is determined 
by the shape of the acoustic wave beam, which depends 
on lotjs, on the aperture and on the focal length of the 
transducer. A focusing ultrasound transducer typically 
yields 1-2 mm lateral resolution in the focal plane when 
used at several MHz. However, localization of the AO ef- 
fect is about 10 times less accurate along the propagation 
direction of the ultrasound. As a consequence, the opti- 
cal and acoustic waves have to be temporally reshaped 
in order to improve the axial resolution. 

Wang and Ku [5] proposed to encode each axial po- 
sition with a frequency-swept acoustic wave but the 
recording time was greater than the speckle decorrcla- 
tion time, excluding any in vivo experiment. As com- 
monly used in standard echography, acoustic bursts are 
a solution to get axial resolution that involves times com- 
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Fig. 1. (Color on-line) Experimental set-up. HWP: half-wave 
plate, PBS: polarizing beam splitter, AOM1, AOM2: acousto-optic 
modulators, FI: Faraday isolator, UT: ultrasound transducer, fc'. 
camera framerate, /j,: laser frequency. In the tank, the scheme axes 
are modified in order to give a view of the inside. 



patible with speckle decorrclation and medical standards 
on ultrasound exposure. The burst technique has bene- 
fited photo-refractive UOT [6,7] and has also been imple- 
mented with digital holographic detection by Atlan et al. 
[8]. Nevertheless, detecting s signals with a camera whose 
minimum acquisition time is in the order of hundreds of 
/is is quite inefficient. Lesaffrc et al. solution [9], known 
as AOCT, is interesting because it returns to continuous 
acoustic and optical waves, performing a good resolu- 
tion though, thanks to a random phase modulation on 
ultrasound and light. Up to now, AOCT technique has 
only been implemented in a photo-refractive holographic 
scheme. In this Letter, we demonstrate the potential of 
AOCT when combined with digital holography. We first 
experimentally study the influence of the random pat- 
tern characteristics on ID AO images. Then the resolu- 



1 



tion of the technique is tested by performing 2D images 
on thick tissue-like phantoms containing small absorbing 
objects. We finally compare two AO profiles of the same 
object respectively obtained with AOCT and the burst 
technique. 

Figure f depicts schematically the experimental set- 
up. A single-frequency 300-mW CW TkSapphire laser 
(Coherent, MBR 110) with a coherence length of 300 
m is tuned at 780 nm in order to work in the optical 
therapeutic window. A polarising beam splitter (PBS) 
divides light into an object beam and a reference beam. 
The object beam, after being amplified by a 2.5 W semi- 
conductor amplifier (MOPA, Sacher Lasertechnik) illu- 
minates a scattering sample immersed in a transparent 
water tank and insonified by a single element acous- 
tic transducer (Panametrics A395S, /[/s=2.3 MHz, fo- 
cal length=78 mm, diameter=38 mm). The transmit- 
ted light passes through a 6-mm diameter circular di- 
aphragm and interferes with the reference beam on a 
12-bit, 1024x1024 pixels fast CMOS camera (Photron 
FastCam SA4) which records images at frame rate fc 
(3.6 kHz maximum frame rate at full resolution). Two 
AO modulators (AOM 1, 2) are placed in each arm of the 
interferometric set-up and shift the frequency of light of 
Iaomi,2- By choosing j 'aom2~ j "aomi = fus + fc/%, we 
exclusively select the tagged photons because the other 
contributions vary too fast to be caught by the camera 
[10] and we perform a two-phases demodulation in order 
to reduce the mean intensity term [11]. The off- axis con- 
figuration adds a spatial separation of modulated and 
non-modulated light in the Fourier domain [12]. 

To get axial resolution in AOCT, the same random 
phase modulation $(i) is applied to the optical refer- 
ence wave and to the acoustic wave via two synchronized 
waveform generators. The reference beam is delayed of 
a time r compared to the acoustic signal in order to se- 
lect a probed zone localized at a given yo position within 
the sample (yo = vus T )- The interference cross term de- 
tected by the camera is the same as in standard AO 
imaging apart from the phase modulation factor. It can 
be written as: Et Er* Gxp[j(j)(t — y / vjj s)] ex P[~ j4>(t ~~ 
t)] +c.c, where Et and Er are the complex amplitudes 
of the tagged photons field and the reference field, vus is 
the velocity of ultrasound (in water vjjs — 1500 m.s -1 ), 
and c.c. is the complex conjugate. Lesaffre et al. [13] 
did the theoretical study of the tagged photons field in 
AOCT in the case of a {0; tt} random phase jump modu- 
lation. They demonstrated that the tagged photons field 
is proportional to the correlation product between the 
two random phase modulations, both in amplitude and 
spatial extent. 

In our configuration of holographic detection, the same 
phase modulation, based on {0;7r} phase jumps occur- 
ring every St, is implemented. The random pattern is 
recorded in the memory of a waveform generator (Tek- 
tronix AFC 3252) which limits the signal length to 524 
fis. As a consequence, the pattern has to be repeated 
several times when the camera acquisition time exceeds 



524 /is. The first experiment aims at studying the depen- 
dence of the AO signal on the modulation characteristic 
time St. The results are compiled in figures 2 and 3. 
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Fig. 2. (Color online) (a) AO profiles along y axis of a scattering 
phantom (thickncss=3 cm, fi' s = 10 cm -1 ), containing a black 
inked cylinder of 10 X 6 X 6 mm? (x, y, z). The same AO profile is 
represented with different resolution values. fc=500 Hz. (b) Shape 
of the absorbing inclusion, extracted from the AO profile with a 
2.25 mm resolution, (c) Numerical simulation of the normalized 
autocorrelation function of the phase modulation sequences used 
to obtain the AO profiles in (a). 

Figure 2(a) shows axial AO profiles obtained in a 
3-cm thick (along z axis) 10%-Intralipid and agar gel 
containing an absorbing cylinder for several resolution 
values. The transport mean free path I* of this phan- 
tom is about 1 mm. The phase modulation pattern has 
an initial length of 200 /is (100 jumps, St = 2 /is) and 
is stretched or compressed by changing the reading pe- 
riod of the waveform generator in order to monitor the 
characteristic time and thus the axial resolution. The 
typical shape of the autocorrelation function of the ran- 
dom phase jump sequence is a triangle of Ay = vus^t 
in width at half maximum, as shown in figure 2(c). The 
theoretical minimum for Ay is one period of the acoustic 
wave, i.e. 0.65 mm at fus=2-3 MHz. When St increases 
from 1.5 to 10 fis, the peak width of the autocorrelation 
function becomes larger so that the absorbing inclusion 
appears in the AO profile with a contrast decreasing from 
95 to 54%. Figure 3 shows the evolution of the ampli- 
tude of the AO signal and the inclusion contrast and 
size as a function of the resolution. These curves are ex- 
tracted from the experimental measurements. In order 
to calculate the characteristics of the inclusion revealed 
by the AO profiles (contrast and size), each profile is 
fitted with a Gaussian envelope used as an approxima- 
tion of the scattering photons intensity profile along the 
acoustic axis. The difference between the envelope and 
the measured profile reveals the inclusion shape as seen 
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in figure 2(b). 
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Fig. 3. Characteristics of the AO profile as functions of the reso- 
lution: (a) amplitude, (b) contrast, (c) full width at half maximum 
of the inclusion (FWHM). 

Figure 3 demonstrates that the AO signal amplitude 
and the contrast of the inclusion vary in opposite direc- 
tions as functions of Ay, which is expected since the peak 
of the autocorrelation function of the random phase pat- 
tern grows in amplitude and widens when Ay increases. 
Thus, the difficulty to perform high contrast AO profiles 
lies in dealing with weak signals. As long as the resolu- 
tion is less than the absorbing object size, the contrast 
and the FWHM of the inclusion are nearly constant and 
give an estimation of the inclusion diameter. If we look 
at the first points on graph 3(c), the extracted diameter 
of the cylinder is 7.5 mm which is 25% as large as the 
real one (6 mm) . This difference can be explained by the 
inhomogeneity of the photons distribution in the pres- 
ence of an absorbing object. Indeed, the probability of 
photon absorption is larger in the neighbourhood of the 
inclusion than in the rest of the scattering medium. 
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Fig. 4. (Color online) (a) AO image of the same phantom 
as in figure 2, performed with Ay = 3 mm axial resolution, 
/c=500 Hz. (b) AO image of a xy section of a scattering phan- 
tom (thickncss=4 cm, // s = 10 cm -1 ) with 2 absorbing inclusions 
of 3 X 3 X 5 mm 3 (x, y, z), performed with Ay = 2.25 mm axial 
resolution, /c=500 Hz. 

An xy AO image of the same sample, presented in 
figure 4(a), is obtained by scanning both the US trans- 
ducer along x and the delay r between the two random 
phase patterns. The quality of this imaging technique 
has also been tested by performing a 2D AO image of a 
4-cm thick 10%-Intralipid and agar scattering phantom 
(I* = 1 mm) containing two black-inked absorbing in- 
clusions separated by 3.5 mm (figure 4(b)). The random 
pattern jump time is St = 1.5 /xs, giving a spatial resolu- 
tion of 2.25 mm. The two inclusions are clearly separated 
and their size is well retrieved. 

The advantage of using random phase jumps instead 
of bursts to achieve axial resolution is the noise reduc- 
tion on the AO signal. In figure 5, the same scattering 
phantom is scanned along the acoustic axis with both 
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Fig. 5. AO profile along the acoustic axis of a 10%-Intralipid 
and agar scattering phantom (thickness=2 cm, fi' s = 10 cm" 1 ) 
with a black-inked absorbing inclusion of 3 X 3 X 6 mm 3 (x,y,z), 
performed with Ay = 2.6 mm axial resolution thanks to random 
phase jumps (filled dots) or bursts (empty dots), /cr=l kHz. 



techniques in the same conditions of acoustic tagging. 
As a consequence, the signal amplitude, the inclusion 
contrast and the FWHM are equal within less than 10% 
on the two profiles. However, the signal-to-noise ratio on 
the AO profile is 3 times worse with the burst technique 
because the useful signal represents only a short frac- 
tion of the acquisition time and is thus surrounded by a 
parasitic signal. 

In conclusion, we have successfully implemented a ran- 
dom phase modulation on US and light to perform dig- 
ital holographic UOT with a millimeter resolution. The 
experiments have been run with continuous acoustic and 
optical waves but the next step is to apply AOCT with 
a 1-ms-pulse laser in order to keep a high level of signal 
while observing medical safety standards. 

The authors gratefully acknowledge the Agence Na- 
tionale de la Recherche for financial support (project 
ICLM-ANR-2011-BS04-017-01). 
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